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ABSTRACT 

Acquiring near-infrared if-band (2.2[im) photometry for RR Lyrae variables in globular clusters and 
nearby galaxies is advantageous since the resulting distances are less impacted by reddening and metallicity. 
However, if-band photometry for RR Lyrae variables in M5, Reticulum, M92, uj Cen, and M15 display 
clustercentric trends. HST ACS data imply that multiple stars in close proximity to RR Lyrae variables 
located near the cluster core, where the stellar density increases markedly, are generally unresolved in ground- 
based images. RR Lyrae variables near the cluster cores appear to suffer from photometric contamination, 
thereby yielding underestimated cluster distances and biased ages. The impact is particularly pernicious 
since the contamination propagates a systematic uncertainty into the distance scale, and hinders the quest 
for precision cosmology. The clustercentric trends are probably unassociated with variations in chemical 
composition since an empirical X-band period- magnitude relation inferred from Araucaria/VLT data for RR 
Lyrae variables in the Sculptor dSph exhibits a negligible metallicity dependence: (0.059±0.095) x \Fe/H\ zw , 
a finding that supports prior observational results. A future multi-epoch high-resolution near-infrared survey, 
analogous to the optical HST ACS Galactic Globular Cluster Survey, may be employed to establish if-band 
photometry for the contaminating stars discussed here. 

Subject headings: globular clusters: general — stars: distances — stars: variables: RR Lyrae — techniques: photometric 



INTRODUCTION 



Photometry of RR Lyrae variables and stars in glob- 
ular clusters is used to establish astronomical distances 
and the age of th e Universe (jSzewczyk et al. 2008: 



Benedict et al.|[2011 ). Reducing uncertainties associated 
with such photometry is a pertinent endeavour in the 
present era of precision cosmology, since even minute 
improvements to the distance scale can facilitate the 
breaking of dege neracies plaguing the selection of a cos- 
mological model (|Macri & Riessll2009h . A mere ~ m .05 
systematic offset in distance for globular clusters arising 
from problematic photometry implies a potential ^5% 
shift in derived ages for the clusters and the Universe. 
Age estimates inferred for globular clusters set a crucial 
lower limit to the age of the Universe, while distances 
for the clusters and their RR Lyrae variable constituents 
are ultimately employed to constrain Hq , from which the 
age of the Universe stems. Fundamental parameters es- 
tablished for RR Lyrae variables are used to anchor 
and calibrate other distance indicators (e.g., the globu- 
lar clu ster luminosity func t ion a n d Cepheid met a llicity 
effect , Di Criscienzo et al. 20061 Maiaess 2010t Feast 
20111 ). which emphasizes the importance of eliminat- 
ing the propagation of systematic uncertainties into the 
cosmic distance ladder. 



Maiaess et al.l ([20121 ) noted that Johnson- Cousins VI C 
photometry for RR Lyrae variables in several globular 
clusters (e.g., M15) exhibit a dependence on the (pro- 
jected) clustercentric distance. The principal (yet not 
the sole) source giving rise to that trend is photometric 
contamination (see also t he discussion in lYannv et al 
1994 IStetson et aT1l2005l ). As the cluster core is ap- 



proached the potential for photometric contamination 
increases in prop ortion to the ste l lar den sity, which rises 
rapidly. Indeed, iMaiaess et al] (|2011bl and references 
therein) asserted that blending/crowding compromises 
the Cepheid distance scale, since variables sampled near 
the cores of galaxies are typically brighter than their 
counterparts occupyi ng the outer regions (e.g., M101, 



Kennicutt et al.ll 19981) . However, the impact of blending 



on extragalactic variable star (i.e., Cepheids) distances 
is disputed, and many researchers instead favour metal- 
licity differences as the origin of the brigh tening (e.g., 
Shappee fc Stanekl l20Llb iGerke et all I2OIII) . Metallic- 



ity likewise increases with decreas ing galactocentric dis- 
tance (e.g., Kudritzki et al. 20121) . thereby introducing 
a degeneracy that needs to be overcome. 

In this study the Maiaess et al.l ( 20121) analysis is 
extended to the near-infrared (~ 2.2 fim) regime, 
and the impact of photometric contamination is as- 
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Fig. 1. — Wesenheit V7 c -band and 7f-band relative distance moduli computed for variables in M15, M5, and NGC2419. The 
moduli exhibit a dependence on position (i.e., declination, S) and clustercentric distance (R), whereby the distances appear 
underestimated for stars near the cluster core. A similar dependence is noted for both jf-band and reddening-free Wesenheit VI C 
distances. That suggests that the former are not being severely impacted by differential reddening. 



sessed for if-band photometry of RR Lyrae v a riable s 



in the glo bular clusters M5 (ICoppola et al.l 120111) 



Reticulum jDall'Ora et ail 12004 iDall'Oral 120051). M9 2 



(iDel Principe et alj|2005l)7a; Cen (jSollima et al.ll2006al) 



and Ml 5 ( Longmore et al. 1990f ). The conclusions pre 



sented do not mitigate the importance of the studies dis- 
cussed, and aspects of the RR Lyrae pulsation analysis 
remain unaffected. For example, the derived pulsation 
periods place seminal constraints on mo dels since pe- 
riod change is a proxy for stellar evolution (jKunder et al 



20111 : Uurcsik et alj|2012l ) 



2. ANALYSIS 

In order to assess the impact of photometric contam- 
ination the residuals of the i-f-band period-magnitude 
relation are evaluated as a function of clustercentric dis- 
tance. The canonical distance modulus for RR Lyrae 
variables is defined as: 

/i = K - M K - A K 
M K = axlogP + f3 

= K - a x logP ~ P ~ A K (1) 

where K is the mean magnitude, Mk is the corre- 
sponding absolute magnitude, logPo is the fundamen- 
tal mode period (overtone pulsators are corrected us- 



ing logP ~ logP + 0.13, ICaputdl2012fl . f3 is arbitrary 
since relative rather than absolute distances are evalu- 
ated to assess the contamination and stars in a given 
globular cluster are assumed to lie at the same dis- 
tance, and Ak accounts for dust extinction in the K- 
ba nd. A period- Wesenh eit (VI C ) analysis was adopted 
by Maiaess et al. ( 20121 ) to negate concerns that differ- 
ential reddening might affect the clustercentric trends. 
The VI C Wesenheit function is comparatively insensi- 
tive to metallicity and also reddening free, but requires 
a photometric colour as input. The infrared photom- 
etry examined here is tied primarily to if-band obser- 
vations, and the effects of differential extinction are not 
addressed directly save by demonstrating commonalities 
with the VI C Wesenheit trends (Figs. [Q [3J. Each RR 
Lyrae variable in a given cluster is therefore assumed 
to suffer the same (unknown) extinction Ak- If siz- 
able differential reddening is present, it likely stems from 
patchy foreground dust distributed non-symmetrically 
across the cluster, and does not invalidate the broader 
conclusions drawn here. Moreover, the impact of dif- 
ferential reddening is minimized in the i f -band since 
A K /E(B-V) - 0.35 (|Cardelli et al.lll989l ). and indeed, 



there are numerou s other adv antages to employing K- 
band photometry (jBonol 120031 ). Nonetheless, in general 
a reddening- free VI C ,JHK Wesenheit analysis is preferred 
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Fig. 2. — Top left, the density of HST ACS stars in M15 found r < 0.8" (well within typical ground-based seeing) from the 
catalogued RR Lyrae variables increases with decreasing clustercentric distance. Top right, the sum of the V-band flux for 
those neighboring stars divided by the mean RR Lyrae variable l/-band flux (Ei**,v/-PiiR,v) likewise increases with decreasing 
clustercentric distance. Bottom, a cropped HST ACS image featuring two RR Lyrae variables near the core of M15. The RR 
Lyrae variables are the brightest stars in the two groups, and the red lines delineate ~ 1". The HST ACS data imply that multiple 
stars lying in close proximity to RR Lyrae variables located near the cluster core are generally unresolved in ground-based images. 
The flux measured for variables in the latter is subsequently contaminated by those stars. 



and desirable to eliminate concerns regarding the impact 
of differential reddening for all clusters. If j3 and Ak are 
chosen to be arbitrary, the above equation can be recast 
as: 



Mo = K - a x log P - @ - A K 
Mo = K - a x log P + 7 



(2) 



The slope (a) adop ted for the period- My function i s ap- 
proximately —2.3 (jCatelan et al.ll2004t lFeastil201lL and 
references therein), and the conclusions derived below 
are not sensitive to that selection (to within the un- 
certainties). The residuals do not display a significant 
dependence on pulsation period. 

J2000 coordinates for the RR Lyra e variables are 



availa ble vi a CDS and the compi lations of lClement et al 
(|200l[) and ISamus et all (|2009l ). For RR Lyrae vari- 



ables discovered in Reticulum, J2000 c oordinates were 
tabul ated using the Aladin en vironment dBonnarel et al, 



2000h and the finder chart of lWalkerl (|l992h . 

In Fig. [1] relative if -band distances calculated for RR 
Lyrae variables in M5 are plotted as a function of po- 



sition (i.e., declination, S) and clustercentric distance 
(i?). Numerous stars in the cl uster appear to suffer 
from contaminated photometry. ICoppola et al.1 (|201ll ) 
remarked that the if-band photometry associated with 
numerous RR Lyrae variables in M5 was likely contami- 
nated by neighboring stars (see also iDi Criscienzo et aL 
20 111 concerning blended objects in NGC2419). The 



extent, scale, and clustercentric dependence of that con- 
tamination is conveyed in Fig [TJ The distance mod- 
uli are systematically underestimated for variables near 
the cluster core, where photometric contamination is 
most pronounced as a result of higher star densities 
(Fig. [2]). PSF photometry does not correct for multi- 
ple stars coincident well within the FWHM of the sys- 
tem. For example, hig her-resolution HST photometry 
(jSaraiedini et al.l 120071 ) demonstrates that the density 
of stars lying within the seeing limit typically associ- 
ated with ground-based images increases with decreasing 
clustercentric distance (Fig. [2]), and thus RR Lyrae vari- 
ables located near the core are susceptible to increased 
contamination by neighboring stars. The contaminat- 
ing flux increases with decreasing clustercentric distance 
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Fig. 3. — AT-band (red dots) and Wesenheit VI C (black dots) 
relative distances for RR Lyrae variables in M15 plotted as 
a function of the clustercentric distance ( R). The if-band 
data were acquired through a large aperture (jLongmore et alJ 
Il990l ). and thus contamination begins to appear at a larger 
clustercentric distance. 



(Fig.©. 

Fig. [U also features the relative distanc e moduli com- 



puted using a VI C Wesenheit function (IMaiaess et al 



20121 ) for RR Lyrae variables discovered in Ml 5 and 



NGC2419. The results inferred from K and VI C photom- 
etry generally follow a similar trend, as demonstrated in 
Figs. Q] and |3l which in part suggests that differential 
reddening is not the cause of the trends noted in the 
-ftT-band photometry. However, it should be noted that 
V7 c -based relative distanc e moduli for RR Ly rae vari- 
ables in M3 are anomalous ( Maiaess et al.|[2012 , and ref- 
erences therein), and display an opposite trend to that 
inf erred for the bulk of the clusters examined here and 
by IMaiaess et al.l ( 2012 ) . Variables near the core of M3 
yield VT c -based distance moduli that are overestimated 
relative to variables near the cluster periphery, for rea- 
sons that are not entirely clear. The acquisition of inde- 
pendent multi-epoch VI C K photometry for that sample 
is needed to help identify the source of the discrepancy. 

Relative distance moduli calculated for RR Lyrae 
variables in Reticulum, M92, and uj Cen are plotted 
as a function of clustercentric distance in Fig. |4j All 
clusters display radial trends whereby variables closest 
to the cluster core are generally brighter than objects 
near the periphery. The distribution of RR Lyrae vari- 
ables in M5, u> Cen, M92, and M15 (K-band) is such 
that a less-biased distance may potentially be estimated 
by relying on stars near the cluster periphery (Figs. [1] 
EJ. However, the trends displayed for Reticulum present 
problems for deriving an unbiased distance. RR Lyrae 
variables in that cluster exhibit small random scatter, 
yet the extent of the pernicious systematic uncertainty 
is unclear (see also M15 Wvi c data, Fig.QJ. Supplemen- 
tal observations are desirable. Nevertheless, the distance 



offset between RR Lyrae variables located near the core 
and periphery may be estimated for M5, M15 (if-band), 
uj Cen, and M92. The resulting estimated offsets are 

Afi 0.09, -0.13, -0.07, -0.06 for the aforementioned 

clusters, respectively. The differences are evaluated for 
illustrative purposes, and are subjective and uncertain 
owing partly to an inability to overcome poor statis- 
tics to assess the extent of the compromised photometry 
(e.g., M92 and M15 (W V i c data)). The bias noted may 
be reduced by remaining cognizant of the trends high- 
lighted in Figs. Q] and [4] and applying robust estimators, 
among which median and sigma-clip algorithms are per- 
haps the simplest, rather than least-squares or average 
which are sensitive to outliers. Indeed, the trends high- 
lighted in Figs. [1] and 0] provide the rationale needed 
when eliminating certain outliers. 

Pearson coefficients computed for M5, M15 (K- 
band), Reticulum, M92, and uj Cen are: r = 0.35 (p 
= 0), r = 0.28 (p = 0.21), r = 0.36 (p = 0.07), r = 
0.37 (p = 0.26), and r = 0.25 (p = 0.03), respectively. 
Pearson coefficients range from r = — 1 to +1 and indi- 
cate the strength of negative and positive correlations, 
respectively, while r = indicates a lack of correlation 
and p represents the significance. The results are to be 
considered in tandem with the evidence hitherto tabu- 
lated, which in sum indicates that the if-band residuals 
for the clusters examined are correlated with the cluster- 
centric distance. However, caution is warranted when 
interpreting the coefficients. Blending is radially de- 
pendent and non-linear, and the blending trends differ 
between clusters owing to the different radial sampling, 
seeing, core and apparent density, cluster distance, etc. 
(Fig. |3|). The coefficients cited are an average for the 
entire cluster, which is typically weighted between the 
marginal to null correlation existing for stars near the 
cluster periphery, and the stronger correlation tied to 
stars near the core (Fig. [1]). 



2.1. 



K-BAND METALLICITY DEPENDENCE 



An argument may be made that the trends displayed 
in Figs. [1] and [4] are tied to a (radial) metallicity gradi- 
ent. Existing estimates concerning the impact of metal- 
licity on if -band distances display a spread (e.g.. iFeast 



201 ll and references therein), although empirical esti- 
mates tend to favour a small dependence. The impact 
of metallicity on distances to RR Lyrae variables as in- 
ferred from if-band photometry is now assessed to com- 
plement e xisting work with an independent determina- 



tion (e.g.. ISollima et al.ll20 06a) 



K s photometry was acquired as part of the Araucaria 
project for 78 RR Lyrae variables in the Sc ulptor dwarf 



spheroidal galaxy ( Pietrzvhski et al. 20081) . The pho- 



tometry can be correlated with metallicities established 
for RR Lyrae variables in t hat galaxy using VLT spec- 



tra (jClementini et al.l [2005) . The metallicity for each 



RR Lyrae variable was determined by comparing the 
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Fig. 4. — Tf-band relative distance moduli for RR Lyrae 
variables in Reticulum, M92, u Cen, and M15 are plotted as 
a function of the clustercentric distance (R). The data are 
mirrored (open circles) to better convey the clustercentric 
trends. The best fit exponential functions are overlayed. 



H and Ca II K line strengths ag ainst a calibration tied 
to variables in globular clusters (jClementini et alj|2005 , 
and references therein). A total of 29 variables were 



found to b e in common between t he IPietrzvnski et al, 
(|2008h and lClementini et al.l (|2005l ) datasets. Such data 
are relevant since the comparatively small stellar den- 
sity of the dSph galaxy reduces biases introduced by 
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Fig. 5. — Ks data for RR Lyrae variables in Sculptor 
(IClementini et al ] |2005l ; IPietrzvnski et al.l [2008l l are plotted 
as a function of the iron abudance. The best fit yields an 
insignificant slope of (0.059 ± 0.095) x [Fe/H] zw ,, and im- 
plies that the K a observations are comparatively insensitive 
to metallicity. The uncertainties may be mitigated by obtain- 
ing K 3 data for additional RR Lyrae varia bles in the Sculp- 
tor d Sph that possess metallicity estimates l|Clementini et al] 
I2005I ). 



crowding/blending. The significant metallicity spread 
also permits an evaluation of the impact of that vari- 
able on distance. The following expression was inferred 
from the data after: the removal of two outliers, con- 
sidering uncertainties in K s , c onsidering a global 0.1 5 



dex uncertainty in [Fe/H] ZH/ ( Clementini et al. 2005), 
and adopting a fixed coefficient for the pulsation term 
(-2.3): 



K s = (0.059 ± 0.095) [Fe/H] 



zw 



2.3 log P + (18.51 ±0.18) 



The results support prior findings that the if-band 
period-magnitude relation for RR Lyrae variables is 
comparatively insensitive to v ariations in chemical 
composition (Fig. [5l see also Sollima et al.l l2006at 



Borissova et al. 200 



Feastll2011 , and references therein). 
Consequently, the trends displayed in Figs. [1] and 0] 
are probably not tied to metallicity. Numerous ad- 
ditional RR Lyrae variables in the Sculptor dSph t hat 
possess metallicity estimates (jClementini et al ] l2005h do 
not have -ftT-band observations. Acquiring near-infrared 
photometry for those stars will reduce uncertainties as- 
sociated with the aforementioned determinations. More- 
over, a holistic approach to constraining the impact of 
metalli city on RR Lyrae and Cepheid distances is pre- 



ferred (|Majaess et al.l l2011bl and references therein) 



rather than relying solely on individual metallicity esti- 
mates for RR Lyrae variables in a given galaxy. A com- 
prehensive analysis is beyond the scope of the present 
study si nce multiple l i nes of evidence sh o uld be em- 
ploye d (jMaiaess et al.l l2011bt iFeasd 1201 ll iMatsunaga 
I2012L and references therein) to assess such an impor- 



5 



tant correlation. 

3. SUMMARY 

if-band (2.2/im) photometry for RR Lyrae variables 
in the globular clusters M5, Reticulum, M92, ui Cen, 
and M15 exhibit a dependence on clustercentric dis- 
tance. Contamination by neighbouring stars causes RR 
Lyrae variables near the cluster centres to appear spu- 
riously brighter, and hence the distance moduli estab- 
lished for the clusters are underestimated (Figs. [TJ [2J 
|4|). That bias in turn introduces a systematic shift in 
the cluster ages, and by consequence on lower-limit es- 
timates for the age of the Universe. The results are 
probably not tied to variations in chemical composition 
among RR Lyrae variables since the if-band period- 
magnitude relation is relatively insensitive to that pa- 
rameter (Fig. [5]), as inferred from RR Lyrae variables 
in the Sculptor dSph. The impact of photometric con- 
tamination may be reduced in part by an awareness 
of the trends exhi b ited in Figs. [H [21 a nd |H (see also 



Yann v et al.l 11994 IStetson et al 



20051: IMaiaess et ail 



20121 ). Blending/crowding may be the principal (yet 



not the sole) source driving the clustercentric trends 
(|Majaess et al.ll2012h . 

Similar results concerning the impact of photometric 
contaminatio n on WV/„ -bas e d RR Lyrae distances were 
described bv lMaiaess et al.l ([20121 ) (Fig. [J). The radial 
trends observed are not unique to a given dataset cited, 
but may be present (to an extent) in most photome- 
try sain£lin^_the cores_| of globular clusters. As echoed 
by IMaiaess et al. ( 20121) . the principal impetus of the 
present research is to increase awareness of contami- 
nated clustercentric photometry in order to reduce the 
propagation of biases into research on RR Lyrae vari- 
ables, especially in the era of precision cosmology where 
identifying and reducing systematic uncertainties is an 
important pursuit. For example, the Reticulum is par- 
ticularly important since it has the potential to yield 
an accurate distance to the LMC from its rich popula- 
tion o f RR Lyrae variables (|Dall'Ora et al.ll2004t IWalkerl 
20111 ). The distance to the LMC remains a cornerstone 



for ef forts to establish precise cosmological parameters 
(e.g., Riess et all 12011 ). Yet establishing unbiased dis- 
tances to RR Lyrae variables in Reticulum appears chal- 
lenging, despite (admirably) small random uncertainties 
associated with the photometric data (Fig. HJ. 

The trends followed by the RR Lyrae variables likely 
promulgate into near-infrared colour-magnitude dia- 
grams a nd hinder efforts to establish proper isochrone 



fits (see Majaess et al. 20121 their Fig. 3). The evolu- 



tionary morphology displayed for cluster stars near the 
core differs from that exhibited by stars near the clus 



ter p eriphery (see also lYannv et al.lll994t IStetson et al 



20051 ) . A degeneracy emerges because the signatures of 
multiple populations, (putative) colour-gradients, and 
differential reddening may be mimicked by photomet- 



ric contamination. As demonstrated here, the acquisi- 
tion of multi-epoch VI C JHK S photometry for cluster RR 
Lyrae variables, which can be used to establish colour 
excesses and reddening-free distance moduli that are 
comp aratively insensitive t o metallicity (see Fig. [51/ N 2.1 [ 
and IMaiaess et al.|[20lTbl 120121 : iFeastJ l201ll and refer- 
ences therein), offers one potential solution for breaking 
that degeneracy. The advent of new mult i-epoch JHK S 



surveys such as the UKIDSS and VVV (IMinniti et al 



20101 : ICatelan et al.ll201ll ISaito et al.ll2012l ). which are 
monitoring globular clusters in both hemispheres , will 
facilitate that effort (see also lMajaess et al.l 120121 ). Yet 
even higher-resolution observations, analogous to data 
from the HST ACS Galactic Globular Cluster Survey, 
may ultimately be needed to provide a comprehensive 
characterization of the near-infrared data discussed here. 
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